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Tlie literature is reviewed concerning the sensitivity of the Dicke radiometer, excluding gain fluc- 
tuations. Discrepancies are pointed out and a new derivation of sensitivity using a Fourier transform 
metliod is used to resolve these discrepancies and to extend the results to radiometers with lossy 
switciies. 

Experimentally it is shown that radiometers using a half-wave square-law, linear-law, intermediate- 
law, or envelope detectors all have a sensitivity equal to the theoretical full-wave square-law detector 
(within the ±20 percent uncertainty of the experiment). 
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1. Introduction 

At the National Bureau of Standards a switching radiometer is used to c()m[)are microwave 
noise sources with a standard liot load [1],^ and a sum-and-difference-correlation radiometer is 
being deveh>ped to compare microwave cryogenic noise sources which uses a switching radiometer 
as a coni|)onent |2|. 

Ahhough tliere is extensive hterature discussing the Dicke radiometer, existing pubhcations 
are difficuh to use. This is partly due to a series of errors and corrections and partly due to 
divergent assumptions, notation, and forms. The literature as a whole has unresolved differences 
between various statements of the sensitivity for the Dicke radiometer. 

The purpose of this paper is to review the literature (sec. 2), derive the sensitivity of the 
Dicke radiometer in a form which includes switch noise (sec. 3), and to compare the derived 
sensitivity with experiment (sec. 4). The sensitivity is rederived, rather than making a further 
correction, for two reasons. First, a statement of sensitivity is needed in a form general enough 
to include the prior literature. This facilitates identifying discrepancies. Second, a basically 
different derivation helps to avoid prior errors and to increase confidence where agreements are 
obtained. 

2. Review of the Literature 

To facihtate comparison between the various authors in the literature, a brief discussion of a 
model switching radiometer is presented. The literature is compared in terms of this model. 
The symbols of this |)aper may be converted to the notation of the other authors by means of 
appendix A. Also contained in appendix A are equivalent statements of predetector bandwidth 
B, and postdetector integration constant r. 



*Sponsored by Advanced Research Projects Agency, Project DEFENDER; ARPA Order No. 515. 
**Microwave Circuit Standards, National Bureau of Standards Laboratories, Boulder, Colo. 80302. 
' Figures in brackets indicate the literature references at the end of this paper. 
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Figure 1. Dicke radiometer. 



2.1. Basic S>vitching Radiometer 

The basic switching radiometer is shown in figure 1. Essentially the power from a noise 
source at the effective temperature T\ is compared with a noise source at effective temperature To. 
It is usually advantageous to operate the radiometer with a known and variable T? which is made 
equal to Ti . The comparison of T\ to To is effected by switching the input of a receiver from 
T\ to To and examining the output of the receiver for a signal correlated with the switching fre- 
quency. Uncorrelated terms (for example, the internal noise of the receiver) do not contribute 
to the dc output of the correlator. 

The noise voltage y\{t) from the noise source at temperature Ti is connected to the receiver 
through a variable coupler of voltage transmission coefficient a\{t). Similarly yo{t) from noise 
source To is fed through a variable coupler with voltage transmission coefficient ao{t). The internal 
noise of the amplifier referred to the input of the amplifier has effective temperature T„. The 
signal Ysit), the input to an ideal noiseless amplifier, is 



ys{t) = a^{t)y^{t)-\-ao(t)yo{t)-\-yn{t) 



(1) 



where y,i{t) is the voltage needed to represent the amplifier noise. Ideally the transmission 
coefficients a] it) and ao(t) vary from zero to one periodically with period Vo^. Typically a]{t) and 
ao{t) are 180° out of phase. The transmission coefficient a]{t) is referred to as the voltage modula- 
tion. The power transmission coefficient p(t) = a^it) is referred to as the power modulation. 
From the modulation ai{t) is generated a periodic reference signal c{t) which then goes to the multi- 
pHer. The reference signal c{t) is referred to as the correlation and usually is sinusoidal or square 
wave. 
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2.2. General Comments 

Eleven papers (as noted in table 1) are liere reviewed which contain stalcnients concerning 
the sensitivity of a Dicke radiomettM-. There is an underlying agreeineiit among eight of these 
and the remaining three retain many of the essential features. In order to recognize this agree- 
ment, various definitions nuist he related (see appendix A) and account taken of the riKxhilation 
and correlation assumed. 

There is an iinpoitant diirerence among the various agiceing papers. The natutc ol this 
difference is essential [o the geneial understanding ol the lileialure, so it will he discussed helore 
commenting on indivichial [)apers. 

The radiometer described by Dicke [3J modulates the noise source by means oi a variable 
attenuator whicli has nearly the same temperature as the noise source of interest. 'I'Ik |)apers 
of (n)ldstein |4, 5], Strom [6, 7], Strum [8], Knight |9], and Johnson [10], collectively referred to as 
(/SSKj, calculate the sensitivity of a somewhat different radiometer. The GSSKJ version of the 
Dicke radiometer modulates the noise source of interest by an ideal noiseless device, which could 
be a variable attenuator at the absolute zero of temperature. This corresponds to T? ^ °K in 
the model introduced in section 2.1. The Dicke radiometer was originally a null device but this 
concept was extended by C/oldstein to include an off-null operation. This introduces many prac- 
tical differences. F^or example, a radiometer operated off-null is not immune^ from receiver gain 
instabilities. The radiometer described by (iSSKJ is a special case of the papers of Kelly, Lyons 
and Root 111, 12|, Colvin [13], and Tiuri [14] collectively referred to as KLRCT. However, GSSKJ 
and KLRGT differ in detail. For example, GSSKJ's results indicate somewhat great(T fluctua- 
tions in the radiometer output due to the input noise source [T]) than the corresponding KLRCT 
results. 

riie origin of the discrepancy lies in the different technique of obtaining the output of the 
square law detector relative to the input. GSSKJ use an autocorrelation method and handle the 
resulting squared terms by means of a theorem usually attributed to Rice [15]. Rice's theorem 
for square law devices is rigorous for stationary (iaussian processes. KLRCT essentially apply 
Rice's theorem directly, neglecting the usually small nonstationary aspect of the signal at the detec- 
tor input. \Ue errors introduced by this simplifying assumption go to zero as the noise sources 
being compared approach the same temperature {T\ = T2), and the receiver is switched between 
the two sources in such a way that the power to the receiver is constant throughout the switching 
cycle. 

The calculations in this paper are as precise as the GSSKJ calculations. Thus detailed 
agreement with KLRCT is expected only when T\ = T2. 

In the last column of table 1, which appears later in this paper, prior calculations of the radi- 
ometer sensitivity are indicated by modulation and correlation type. Those calculations which 
differ with this paper are marked with an asterisk. Although (^oldstein [4, 5], Strum [8], and 
Knight [9] are marked with an asterisk, they are considered to be in the group of eight having under- 
lying agreement. More will be said of this in the detailed comments which follow. 

2.3. Specific Comments 

In 1946, Dicke [3] suggested modulating a noise signal under study as a means of distinguishing 
the signal of interest from the noise originating in the amplifying system and also to reduce the 
fluctuations in the radiometer output due to receiver gain instabilities. He derived the sensitivity 
of a radiometer using a statistical approach for square-wave modulation and sinusoidal correlation. 
In the analysis Dicke assumed that the detector was square-law, then he generahzed to the linear 
detector. Selove [16] corrects an error in the generalization to the linear detector and Strom 
[6, pp 37 and 68] challenges an overall factor. For a square-wave modulation and sinusoidal 
correlation, Dicke obtained a sensitivity which is a factor of 2 greater than the result to be obtained 
in this paper and thus comparable with the accepted sensitivity of an ideal (no gain fluctuations) 
total power radiometer [11, 12, 13, or 14]. Because the switching radiometer cannot have com- 
parable sensitivity, Dicke's result must be in error. 

129 



In 1955, Goldstein [4] analyzed the switching radiometer using sinusoidal-voltage modulation 
(or switching), sinusoidal correlation and square-law detection making use of the autocorrelation 
function and assuming the noise signal had a (iaussian distribution. However, as noted al)()ve, 
he assumed a noiseless modulator which corresponds to assuming 7^2 = °K. Goldstein's result 
as corrected by himself [5] and Johnson [10] agrees with a special case of the analysis to be 
presented. 

Also in 1955 Bunkin and Karlov [17] analyzed the switching radiometer using the autocorrela- 
tion method. Their analysis assumes a Gaussian noise distribution, square-wave modulation and 
sinusoidal correlation. In addition, the effects of amplifier gain fluctuations are considered. 
This paper depends to some extent on other papers in the Russian literature which are not readily 
available. It was not estabhshed for certain if the paper agrees in all details, but no obvious 
discrepancies occur. 

Strom [6, 7] shows that the sensitivity of the switching radiometer is unaltered if a biased diode 
is used as the detector instead of a square-law detector. He obtains a sensitivity that is less by a 
factor V 2 than the expected value because the out-of-phase rejection of noise by the correlator 
was not taken into account. 

Strum [8] discusses the effect of fluctuations in various components of a radiometer system 
and discusses various types of linear detectors. He comments on the use of an IF amplifier in an 
analysis. Strum defines a generalized bandwidth which reduces to the usual one for rectangular 
bandpass but is not equivalent in general to that used by Dicke: Kelly, Lyons, and Root; Ward and 
Richardson; Tiuri; and in this paper. 

Kelly, Lyons, and Root [11, 12] analyzed a general comparison radiometer (T^i and Tz arbitrary) 
for arbitrary modulation and correlation. They treated gain fluctuations and optimized the modu- 
lation wave form depending on the rms amplitude of the gain fluctuations. They also presented 
a proof that a square-law detector is superior to any other detector law for Gaussian noise (ignoring 
gain fluctuation of the amplifier). Except for the differences off-null already noted, there is no 
disagreement with this paper. 

Colvin [13] discusses the switching radiometer under rather general conditions as part of an 
extensive survey of various radiometers. Again, except for the off-null differences discussed, there 
is no disagreement with this paper. 

Knight [9] considers a radiometer with asymmetrical rectangular modulation and correlation. 
Tlie calculation is of the Goldstein type. In his paper, eqs (3) and (4) lack the factor 27r on the 
right-hand side which occurred in his prior equations. Knight's corrected equations are a special 
case of the analysis to be presented. However, his extension to the balanced radiometer operation 
did not take into account the fluctuation of the output due to the noise sources. 

Ward and Richardson [18] consider a power sinusoidal modulated radiometer with IF ampli- 
fiers of different image and signal responses. The sensitivity obtained is that expected for a total 
power radiometer, and for the same reason as noted for Dicke's paper, this must be in error. 

Johnson [10] corrected (/oldstein's results and extended the calculation to include various 
combinations of sinusoidal and square-wave modulations and correlations. However, in Johnson's 
extension to square-wave modulation his (1/4) crlcr^ term should be replaced by erf erf,. ^ 

Tiuri [14] states the radiometer sensitivity for several modulations and correlations. The 
results are the expected ones within the approximations already discussed. 

3. Analysis 

The sensitivity of the switching radiometer depicted in figure 1 is determined by calculating 
the ratio of the dc output of the radiometer y/. 2 and the rms fluctuations at the output, ct/. 2. The 
time dependent noise voltages will be denoted by y{t) with a subscript denoting the point in the 
circuit shown in the circles in figure 1. An W amplifier is not used in figure 1 but is treated in 
appendix C. 



^Tliis lias Ix't'ti acknowlcdjied l)y Johnson in a private communication. 
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3.1. The Fourier Transform 

The Fourier transform preserves required phase iiiiorniation. Because the transform of a 
random function does not converge, all the noise voha«j^es y(t) are defined to he zero except for 
— TI2<t<TI2. The liuncatt^d function of y{t) has identical [)r(>perties with y{t) in tlie hmit 
T^^. Thus by y{t) it will he understood the truncated function oi y{t) . The Fourier transform 
of y(^) is denoted F(f), again with a subscript to denote [)oint in the circuit. 

The foUowing form is used for the Fourier transform |)air: 



(2) 



y{t)= I " eJ^'""Y(v)di'. (3) 

J — oc 

The Fourier transform or inverse transform of a product (in either t or v space) is the convolution 
of the transforms of the two factors (convolution theorem). With the choice in (2) and (3), the 
convolutions (denoted by it) are defined by 



yAt)iryAt)^\ yAt-t')y,[t')dt' , (4) 

and 

YAv)-kYAp)^\^ YAv-v')YAv)dv'. (5) 

3.2. The Power Spectral Density and the Wiener-Khintchine Theorem 

The power spectral density W [v] is related to the Fourier transform oi y{t) in the following 
manner [19], 

JrM = T-^oo Y , (6) 

where the asterisk denotes the complex conjugate. It can be shown that W{v)dv is the i)ower 
dissipated in a 2-fl resistor due to the spectral components between v and v-^ dv in the voltage y[t). 
The Wiener-Khintchine Theorem is stated [20] 

W{v)={^ R{t)e-J^^'fdt, (7) 

so that 

R{t)=i^ W[v)e^^^^'dv, (8) 



where R(t) is the autocorrelation function defined by 



RiD^V-^ ^ {'^'yinyit'^Ddt', 

^ J-T/2 



(9) 



Using (8) and (9) for / = 0, 



yHt)= W{v)dv, (10) 



where the bar indicates the time average in the interval —7/2 to + r/2. 
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The analysis is ultimately stated in terms of measurable properties of the input noise signals 
[21,22]: 

WM=lkTu 

W2{v)=\kT,, (11) 

Wn{v)=\kTn. 

where k= L3805 X 10^^^ J °K"^ is the Boltzmann constant. The factor 1/2 in (11) occurs because 
the power spectral density is defined for both positive and negative frequencies. 

3.3. The Switch 

A model for the switch discussed in section 2.1 led to eq (1). For this analysis a more general 
y^t) composed of an arbitrary number of periodically modulated noise signals will be used. The 
voltage modulation of the ith signal [a\{t)] is analyzed into its Fourier components: 



ai{t) = ^Aj,n,e^^^''^'^\ 



(12) 



where m takes all integer values between —^ and +00. Similarly the correlation signal is 
decomposed: 

c(0 = 2^'"^''"''''"'- (1^) 

The following associated transmission coefficients are also required: 

Pi[t)^d][t), (14) 

qi{t)^p]{t) = att), (15) 

Uij(t) = ai[t)aj{t), (16) 

Vij{t)^af{t)aJ{t), (17) 



which are decomposed into 



and 



m 

9'U) = 2 <?'■-'« «''"'""'°'' (19) 

m 

uij(t) = ^U,j,,ue^"'^^'"'*, (20) 

m 

i;y(0 = 2Fij,„,e-''«2w_ (21) 

m 
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3.4. The y(t), Y(p), W(v) Description 

Each element of the radiometer depicted in fi<2;ure 1 is described hy the modification to y(t) 
before and after the element. Thus, 

ys{t) = ^aM)rAt) (22) 

X4(0= /i.i(0*r.v(0, (^'^) 

yiAt) = y\{t\ (24) 

y.l(0=/^^lW*rD(0, (25) 

yM{t) = c(t)yyx{t), (26) 

yF2{t) = h,.{t)i^yM{t). (27) 

where hA[t), hy^{t), hf^it) are the impulse response functions [20J ()f the amplifier, bandpass fiher 
and low [)ass fiher, respectively, and c{t) is defined in (13). 

Taking the Fourier transforms of (22) to (27) and using (12), 

YsM=^ Ai, -,,F/(^ + m^o), (28) 

Y.-M = H.MYs(i^), (29) 

Y„{i') = YAt')ifYAt'), ... (30) 

Y,.;{v) = H,,{i')Yi>{i'), (31) 

Yxi{v)= J C-,„Yr,{i> + mp„}, (32) 

//J = - 00 

Y,'2M = HMi^)Ym{p). (33) 



The Fourier transform of the impulse response function is known as the system function, or 
transfer function, denoted H{p), 



(.)-/_ 



H{v)^ e-'''^''h{t)dt, (34) 



and will be used to define the power response function Gip) 



G{v)^H[v)H''{v). (35) 
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The power spectral density is obtained from (28) to (33) by means of (6) and (B-1) of appendix 
B. This results in 

WAi') = Ga{i')Ws{v), (37) 

Wo(v)=\'^^T-'Yo(i^)YUi^), (38) 

Wfiiv) = GFi(v)Wo{v), (39) 

W.u{v)= J C-„,C*,„Wri{i'+mvo), (40) 

m = - X 

Wt2{v) = GF2[v)WM{v\ (41) 

where Ga(v) (see eq 35) is the power response of the amplifier, etc. The spectral density of Woiv) 
will not be deduced from Wa{v) due to the complication arising from Wa{v) being nonstationary. 

3.5. Calculation of yy-i 

The calculation of yF2 (d-c output of the radiometer) is based on the following relationship 
(see appendix B): 

Lim F, 2(0) 

yF2 = T-.oo J ' (42) 

Using (31) to (33), 

YM0)=Hr2{0) J C-nMFi{mPo)Yo{mvo). (43) 

m=-oo 

The task now is to obtain Yoimvo) in terms of the known power spectral densities of (11). The 
approach is to cast YoimPo) in a form to utilize the relationship (see appendix B) 

Lim F..(0) _ 

r->oc — j; — Aj-, (44) 

where 

K-^j^ WAv)dv, (45) 

and Fj.x(0) is the Fourier transform of y|(0 for ^' = 0, where yx{t) is a noise function with zero 
mean. 

If we assume a flat amphfier response so that Ha{v-\- mvo) —Ha{v), then (29) may be written, 

Ya{v) - 2^'-' -pY^i.A{v^pv^), (46) 

up 

where Y\ifA{v) is the Fourier transform of /z,i(^)^yi (0 , etc. The approximation for the transfer 
function of the amplifier, Ha{v) might seem unrealistic for large p, but generally Ai^p decreases 
rapidly with increasing |p|. 
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Substituting (46) into (30) will <2;ive terms of the form 

y,(^ + />^o)•F,(l^ + (7^^o) = yr^(^+/^A^ + ^/^^o), (47) 

where Yxz(i^) is the Koiiiier transform of yx{t)yz{t) , x, z are any of the siil)s(*ri[)ts \^A, 2^ A or 
ni^A. Equation (47) is valid even when;i: = z. For a proof of (47) see appendix B. 

We are now in a position to use (44) to facihtate the evaluation of (30). Because Yj-zi^-^- pPo 
-\-qvo) gives no contribution except for x = z and v^pvo-\- qPo = (see appendix B), then by re- 
taining only nonzero terms arui excluding m = (because Aj, o = 0) , 

Urn Ypimvo) _^ 

T^^ j; —2^Ai,qAi,-q-rnKii,A, (48) 

i, q 

where from (6), (11), and (45), for i=l or 2, 

K,i,A = I kTi j^^GA{p)dp. (49) 

Thus, for the usual case that Co = 0, (42), (43), (48) and (49) combine to give. 






(50) 



where 



and 



X = | m,2(0)[2G,(0)*G4(0)]^/2, (51) 



B 



I GA{v)dv 



2G..i(0)*a4(0) ' 



(52) 



or using appendix A, 



\^Ga{v 



)i^GA{v)dv 
^^ 2G^(0)*G^(0) ' ^^^^ 



and Pi, m defined in (18) is introduced by means of (B— 7) of appendix B. The term B is the con- 
volution bandwidth of the amplifier. For a rectangular power response, B is the width of the 
nonzero response in hertz (cycles per second). 

3.6. Calculation of o"/ 2 

The rms fluctuation (standard deviation) 0-/2 is defined by 

2 

The bar indicates the time average in the interval —7/2 to Tj2. Using (10) this may be written, 

o"?2= WF2(v)dv — yF2. 
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If the 8-function term [23] from (38) is separated from Wf2(i') of (41), denoted by a prime, 

o-F2=f" W',^(v)dv. (54) 

J —00 

Combining (22), (23), and (24), and assuming the modulation frequency is small compared 
with the amplifier bandwidth so the modulation factor distortion is negligible, then 



yz)(0 = 



^ai{t)yii,A{t) 



where 

ya^Ait) ^hA{t)-kyi{t). 
Expanding (55) and using (14) and (16), 

ynit) =^Pi{t)y\i,A{t) ^2^u;j{t)y,^A[t)y,^A{t) 



(55) 



(56) 



(57) 



Taking the Fourier transform of (57), then using (6) and (B-2), to obtain the power spectral density, 
W'^(v) = 2 Pi, -mPl _,„r',^,, ,^,(v + mvo) + 4 2 U;i, -,nUl _,„W\^,, ,*, {v + mvo), (58) 

/, m i<j, m 

where WurAjirA (^) is the power spectral density of r.*^(0y;*.4(0 and the primes indicate 8-func- 
tion terms are omitted. It can be shown [15] that 



where Ga{i^) is introduced in (37). Also for i y^j, 

mi^A, i*A (»') = \_Ga{v)WAv)\ • \_G,{v)W-^y)\ . 



(59) 



(60) 



Using (59), (60), (10), (B-8), and (B-9) to simplify (58) and then using (39), 

Gm(j')[G,(i^)*G,(j^)], 



W',lv) = -k^ 



^Qu„Tf+2^Vu,oTiTj 



(61) 



where Qi, o and F/;, o are defined via (19) and (21). 

For the computation of Wf2{v) using (41), only Wm(v) for v — Q will contribute because of the 
narrow low pass filter F\. Consequently in using (40) and (61), GF\{v±mvQ)[GA{v±TnvQ)^GA{v±mvQ)\ 
may be replaced by GF\{mV{^[GA{raV{^^GA{rnvQ)\ because the power transfer functions are even in v. 
Furthermore, in the case where Cn decreases rapidly with increasing |n|, then to a good approxi- 
mation GA{nivo)i^GA{mvo) — GA{0)i^G{0). Therefore, 



l^; (v) - I XW2 



^QLoTf+2 2Vu,oT>Tj 



GrzMIGr2(0), 



where X is defined in (51) and 



M = 



2 C^mC-mGn{mVn) 

HI 
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1/2 



(62) 



(63) 



Thus using (54), 



^here 



fr,->-XMT '/2 



Xa-,orf+2;^F,-.xT,- 



1/2 



< 



Gr-Ai^)ihlG,-A()). 



(64) 



(65) 



3.7. Sensitivity 

A measure of a radiometer's sensitivity is (r^.^{dyj:^ldTi). It is customary to specify instead, 
tlic minimum resolvable change in Ti, denoted AT^in- Using the criterion that a change in yp2, 
denoted Ayfi due to a change in Ti is resolved when Ay^i ^ o-ri, then 



Using (.50) and (64), 



""" dyF2ldT, 



(66) 



M 



AT„ 



^Q,nrf + 2^V,,,J)Tj 



1/2 



VB^{2C-,„/y,,(mr„)P,,-,„}' 



(67) 



where 5, r, C-^„, ///<,(i^), (?,, 0,^0,0, M are defined in (53), (65), (13), (34), (19), (21), and (63) respec- 
tively. 

a. Radiometers With Premultiplier Filter 

Most radiometers have a narrow (compared with p^) bandpass fiher (denoted Fl in figure 1) 
preceding the muhipher. For this type of radiometer the sensitivity is independent of the cor- 
relation wave form. This follows from the fact that Hi.\{mP{)) in (50) is zero except for m = ±l. 
For this case, 



AT„ 



i i<j : 

V2B7|P,,i|cosi// 



(68) 



/here 



The sensitivity for particular modulation wave forms are shown in table 1. 

b. Radiometers With No Premultiplier Filter 

To simplify the calculations, assume the transfer function of the bandpass filter is constant 
in phase and amplitude and equal to /// i(^'o). In this case (67) becomes 



•2(f^^l|2 



[c'm 



Ar„ 






1/2 



[c{t)pMWB7 



(69) 



where c(t)p]{t) is introduced via (B-10) of appendix B and c^(0 is introduced in an analogous manner. 
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TabLE 1. AT,„i)ifor special ideal switching radiometers. The amplifier noise is 7^3; ^3, o=l for all radiometers; d is the 
fraction of the time the transmission coefficient differs from zero. 



1/ 



Modulation 


Correlation 


F 


Quo 


<?2,« 


Vn, 


r,3.0 


^23,0 


Prior calculations 


1. Square wave 

2. Square wave 


sinusoidal 
square wave 


77/V2 
2 


1/2 
1/2 


1/2 
1/2 






1/2 
1/2 


1/2 
1/2 


a*ce*tg"j*''"t k" 
e*tf>h*1:jni:k" 


3. Voltage sinusoidal 

4. Voltage sinusoidal 


sinusoidal 
square wave 


2 V2 


35/128 
35/128 


35/128 
35/128 


3/128 
3/128 


3/8 
3/8 


^3/8 
3/8 


b*t1:d*t1:e*tjn 


5. Power sinusoidal 

6. Power sinusoidal 


sinusoidal 
square wave 


2 V2 

77 


3/8 
3/8 


3/8 
3/8 


1/8 
1/8 


1/2 
1/2 


1/2 
1/2 


g() i* 1,0 


7. Power sawtooth 

8. Power sawtooth 


sinusoidal 
square wave 


V27r 
4 


1/3 
1/3 


1/3 

1/3 


1/6 
1/6 


1/2 
1/2 


1/2 

1/2 


gf 


9. Rectangular , 

10. Rect.(^^^) 

11. Rect.(^^l) 

12. Rectangular 


sinusoidal 
square wave 
square wave 
rectangular 


7r/(V2 sin ird) 
d-' 

[d{i-d)]-^i^ 


d 
d 
d 
d 


l-d 
1-d 
l-d 
l-d 








d 
d 
d 
d 


l-d 
l-d 
l-d 
l-d 


hn 



*Results of author conflict with this paper as noted in footnote. 
°Resuhs of author agrees in detail with this paper only when T\ =Ti. 
tThe calculation made was equivalent to assuming 72 = °K. 
XA. rectangular filter bandpass assumed in calculation. 

^ Dicke [3] obtains a sensitivity greater by a factor 2 (thus Ar^in less by factor 2) than this paper. 
^ Goldstein [4, 5], a minor algebraic error corrected in ref. j. 
^ Bunkin, Karlov [17]. 

^ Strom [6, 7] obtains a sensitivity less by a factor V 2 than this paper. 

e Strum [8]. The bandpass defined does not agree with this paper except for rectangular bandpass. 
'Kelly, Lyons, Root [11, 12] agree only when radiometer balanced. 

** Colvin [13] agrees only for balanced, constant power radiometers. ' 

^ Knight [9]. Equations (3) and (4) of Knight's paper omit a factor 27r. In addition. Knight's results differ from this paper in that the terms in T\ and T2 do 
not occur. 

' Ward and Richardson [18] obtained a sensitivity greater by a factor 2 V2 than this paper. 

^ Johnson [10] algebraic error for square wave modulation. 

''Tiuri[14]. 



4. Special Results 

A variety of particular radiometers have been considered in the Hterature (see sec. 2). These 
radiometers were hsted in table 1. The radiometer types shown in table 1 are ideal in three ways. 
First, the switching element introduces no noise of its own, second, the maximum transmission 
coefficient is unity, and last, the relative phases of switching and correlation wave forms a{t)^ b{t), 
and c{t) are selected to minimize ^Tj^m- The factor F is equal to 2~*/^|Pi, i|~^ for sinusoidal cor- 
relation or for radiometers with a narrow band precorrelator filter (and any reasonable correlation). 
These radiometers are denoted as sinusoidal correlation in table 1. For other radiometers, F is 
equal to [c^{t)]^l^[c{t)pi{t)y^. Explicitly the power modulation wave forms used are 

(i) square wave modulation, 



P.W= 2 



mn 



. ^m2nvot. 



(70) 



(ii) voltage sinusoidal modulation, 



p\{t)= 3/8 + 1/2 cos 27r vot + 1/8 cos 47r vot; 
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(71) 



(iii) power sinusoidal modulation, 

pi(0=l/2-hl/2cos27r ^^0^; (72) 

(iv) power sawtooth modulation, 

p, (0 = 1/2 +7 2 (277^)-'e'-*-2"'»'; (731 

(v) rectant2;ular modulation, 

^^^^^_ sin^^ (74) 

where ^ is the fraction of the time that pi(t) differs from zero. Except for voltage sinusoidal 
modulation, p2(0 is chosen such that pi(^) +p2(0 =^ 1- 

4.1. Nonideal Switch 

The theory developed applies to arbitrary switching wave forms. There are two common 
deviations from the ideal forms used for table 1. First, the nonideal power transmission coefficient, 
p'iit) of the switch varies from zero to p/, max instead of from zero to one. Thus, 

P/(0=P/,maxP/W, 

so that 

Qi,0=Plm-aJ?i,ih '(75) 

and 

V'ij^ 0= Pi, maxPj, max^u, 0. (76) 

Equations (75) and (76) may be used in place of the corresponding ideal parameters, to extend the 
results of table 1 to the switch with nonideal transmission coefficient. 

A second complication is that the nonideal switch radiates thermal energy in direct proportion 
to its temperature, and the fraction of incident energy absorbed. For a matched three-port device 
at temperature T^^, the effective temperature of the radiation Psw(^)^sw is [24] 

Psw(0nw=[l-PiW-p2(0]nw, (77) 

where the piit) are given by (14). The thermal noise from the switch may be treated formally as a 
noise source of the temperature of the switch which is power modulated by amplitude p^^it). Thus 
for the ideal square-wave modulation, rectangular modulation, power-sinusoidal modulation, and 
power-sawtooth modulation of table 1, 

Psw(0=0. (78) 

For ideal voltage sinusoidal modulation, 

Psw(0=^-^COS ilTVot. (79) 
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Table 2. Balanced, constant switched power, nonideal radiometer 
sensitivities. 
The maximum power transmission coefficient through the switch is 
Pi, max — P2, max — *^; T\ — Ti\ the effective amplifier noise is Tn\ 
and the physical temperature of the switch is T^y^,. 



^T^ln = F 



ri + g 'n+(a^-i)nv 



A. Narrow band precorrelator filter or sinusoidal correlation 



Modulation 



Square wave 
Power sinusoidal 
Power sawtooth 
Rectangular* 



77/V2 = 2.22 
2V2 = 2.83 

V27r = 4.44 
7r/(V2 sin rrd) 



B. No precorrelator filter 



Modulation 


Correlation 


F 


Square wave 
Power sinusoidal 
Power sawtooth 
Rectangular* 
Rectangular* id ^ h) 
Rectangular* {d ^ h 


square wave 
square wave 
square wave 
rectangular 
square 
square 


2 

7T 

4 

[d{l-d)]-''-^ 

d-' 

a-d)-' 



*d is the fraction of the time the transmission coefficient is dif- 
ferent than zero. 



4.2. Balanced Radiometers With Nonideal Switches 

The expression for ATj^in simplifies when the noise sources being compared have nearly the 
same effective temperature and when the power to the amplifier is constant [i.e., pi(t) -{-pzit) 
= constant]. Table 2 lists the results for radiometers having nonideal switches that satisfy these 
special conditions. The sensitivity indicated is for the correlation phase which minimizes Ar^jn. 

The Ar^in for the voltage sinusoidally modulated radiometer does not simplify for Ti = Ti. 
However, for a radiometer with a narrow band precorrelator filter (or sinusoidal correlation) and 
a switch where the maximum power transmission coefficients /?i, max—P2, max — ct^ ^T^^^ for Ti = T2 
is bounded by the following expressions: 



Ar^in > 2 V2 



Tr. + a-Ts+ia-i-^irsw 



(80) 



and 



Ar^jn < 2 \/2 



-r, + a-T3 + (a-i-7)rsw 



(81) 



For no precorrelator filter and square wave correlation, 2 V2 is replaced by tt. 
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5. Experimental Procedure 

The sensitivity of a Dicke radiometer similar to W(^lls, Daywitt, and Miller [IJ was deter- 
mined experimentally. A switchable circulator was used for the switch instead of a magic tee 
and motor driven attenuators. The switchable circulator provided square-wave modulation at 
1000 Hz. The 30 MHz output of the IF amplifier was attetuiated by various amounts and then fed 
to a crystal detector. The detector law of this detector is shown in h<;ure 2. V\w arrows indicate 
operating points for which the radiometer sensitivity was measured. The numbers near the 
arrows are the minimum resolvable temperature changes measured for each level. 

The mixer-preamplifier-IF amplifier voltage response is shown for one chaiuicl in figure 3. 
The other channel is similar. The markers along the bottom of the figure are synthesized at 2 
MHz intervals. The convolution bandwidth was calculated using a computer program based on 
(53). The system bandwidth is determined from a single channel response as discussed in ap- 
pendix C. 

The time constant for the low pass filter of the lock-in amplifier nominally set at 10 sec inte- 
grating time is based on the voltage response measurement shown in figure 4. 

The recorded data were plotted on probability paper. The ordinate of this paper is the am- 
plitude (voltage) of the radiometer output, and the abscissa is the percent of the time the recorder 
trace exceeds the ordinate value. The scale of the paper makes the plot a straight line if the 
data have normal (Gaussian) distribution. From the plot, the mean of the data is the value of the 
ordinate for which the abscissa is 50 percent. The standard deviation a (for normal distribution) 
is the difference in the ordinate values for abscissa values of 50 percent and 84.1 perceht (or 15.9% 
and 50%). Figure 5 shows two sets of data, one before and one after a change of 0.0087 dB (21.6 
°K change from 10,800 °K) of input power is made in one noise source. 

The various parameters for the radiometer under test are given in table 3. The theoretical 
^^min is based on table 2 for the balanced radiometer using square-wave modulation. The 
experimental ^Tmin is the mean result of AJ'i/(Ay/.-2/cr/,2) taken from four plots similar to figure 5 
for a Ar,-2L6°K. 

A feature common to all data and shown in figure 5 is an increase in (j when the noise power 
decreases from null. No change in cr is noticeable when the noise power increases (out to 0.2 
dB change). The cause of this feature is not known. 

The radiometer sensitivity as noted in figure 2 for the radiometer at various operating points 
are comparable. 

6. Conclusions 

The areas of agreement within the literature have been ferreted out and the nature of dis- 
crepancies determined. A new analysis of the Dicke radiometer was presented which was general 
enough to include all of the papers reviewed as a special case and to extend the sensitivity calcu- 
lations to include radiometers with a lossy switch. 

The most sensitive radiometer is a square-wave modulated radiometer with wide band (i.e., 
no premultiplier filter) square-wave correlator. However, the potential 11 percent advantage in 
sensitivity for wide band square-wave correlation over systems using a narrow premultiplier filter 
must be weighed against the practical difficulties associated with wide band systems. Also the 
narrow band correlator insures that the output of the radiometer is proportional to the difference 
of input powers [i.e., to Ti — T>; see eq (50)] independent of the asymmetry of the modulation. 

The sensitivity of a square-wave modulated radiometer using a narrow precorrelator filter was 
determined experimentally. The measured sensitivity was the same as the theoretical radiometer 
using a full-wave, square-law detector. The sensitivity of the radiometer was independent (within 
the experimental error) of whether a half-wave square-law, linear, or intermediate-law detector 
was used. 
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Figure 5. Data plotted on probability paper. 



Figure 3. Amplifier voltage response. 



Table 3. Parameters of the Dicke radiometer under test 



Predetection bandwidth (B) 

Postdetection integrating time (r) 

Receiver noise temperature {T„) 

Maximum switch transmission (a) 

Switch temperature {Tsw) 

Noise sources {T\ — T-i) 

Switching frequency (z^o) 

Switching waveform 

Microwave frequency (local oscillator) 

A^min (theory) 

^^min (experimental) 



9.2 MHz ±10% 

18 sec ± 10% 
750°K±5% 

(2.57)-i±4% 

295 °K 

10,800 ±2.5% 
1000 Hz 
Square wave 
9 GHz 

2.3 °/C ± 15% 
2.3°K±20% 
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7. Appendix A. Equivalent Statements for B and r and Notation Conversion 

for Literature 

Various statements are used in the literature for B and r. This section gives the deli nil ion 
used in this paper in alternate forms which are easier to compare with the literature. 

Using the Wiener-Kliintchine Theorem, (7), and assuming the order of integration can be 
interchanged, the following relationships can be proven: 

C(0)*G(0) = p OMi^, (A-1) 

p G{v)-kG{v)dv=\{^ G(v)dv \ (A-2) 

g(0) - r G{v)dv, (A-3) 

J -00 

//(0)=J h{i)dt, (A-4) 

and Parsevals relation, 

p H[v)n-^[v)dv=[^ h^{t)dt, (A-5) 

where g{t) is the Fourier transform (autocorrelation function) of the power response G(v). and 
h{t) and H{v) are a Fourier transform pair, with h{t) a real function. 

If instead of using frequency v in hertz (cycles per second) as a parameter, the angular fre- 
quency 0} = 27TV is used, then the following conversions may be used: 

v = o)I{2tt), 

Gi,(p) = 27rGco{o^), 

WM = 2ttW,M. 

n(r) = 27rFjco), 

//,(^) = 27r//jco), 

GAv)iKGi,{v) = 2TTGu,[(j})i^Gjcx}), 

where the subscripts v and co are introduced to distinguish a difference in functional form. This 
subscript will now be dropped. The functional form appropriate to the function used is under- 
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stood. The conversion from W{p) to JVioj) is chosen so that 

W{v)dv = JF{a))d(xj. 
In the (o notation, 



and 



F(a,)= (277) 



/: 



e-J''^y{t)dt, 



r(w)=}^i'2 (277/7)yMy*M. 



The conversion system above is not unique. The easiest way to determine what system a given 
author is using is by examining the form used for the Wiener-Khintchine theorem, or statements 
which include the autocorrelation function. The autocorrelation function is universally defined 
in a manner equivalent to (9). 



7.1. Equivalent Statements for B 



The definition of B is in (52): 



B = 



p GAiv)dv 



Using (A-2), 



B = 



/: 



2G.4(0)*G.4(0) 



GA{i>)ifGA{i')dv 



(A-6) 



or (A-1), 



B- 



2G.4(0)*G^(0) 



p G,(i^)c/J' 



(A-7) 



J — X 



(A-8) 



Beginning with (A-8), noting that Ga{v) = Ha(v)H^{p), and converting to w notation, 



B=^-^ 



Aiw)\^da) 



4-7T 



r \ham\ 



(A-9) 



*d0) 



Beginning with (A-8), applying (A-3) to the numerator and (A-5) to the denominator, 

«-'^(0) 



B- 



'j'j'' 



(A-10) 



2 I g'j,U)dt 
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Using a normalized autocorrelation function, pA{t) = gAit)/g.\{0), 

J - oc 

7.2. Equivalent Statements for r 

The definition of r is from (65): 



(A-11) 



-/; 



Gi.-2{v)di'IG,.-,{()). 



(A-12) 



Noting that Gf2{^) = ///«'2(^')/^*2(^) ^md converting to the oj convention, 



|///.2(0j)|2rf0J 



277//2,(a; = 0) 
Using (A— 4), and (A-5), and assuming a causal system such that h{t) =0 for ^ < 0, 



(A-13) 



/: 



hl.it )dt 



I 



hy.iOdt 



(A-14) 



7.3. Literature Conversion 

Several representative statements from the hterature referring to the minimum resolvable 
change in temperature for a radiometer near balance (Ti = To) and for Ar^in ^ Tn will be converted 
to the notation of this paper. 

a. Dicke [1946] 
Dicke's [3] equation (21) for square-wave modulation and cosine correlation is, 



AT^TT 

T 2 



N 



foe 1 r r°° 11 ^/2 

\F{o))\\io} |S(co)|2r/co 



'/ 



(A-15) 



S(0) \F{M)\Hoi 



Using (A-9) and (A-13) and the following conversions, 

T=Tx = T2. 

]\ = {T,-^T,)IT, (for r, = 290 °K), 

F{C0) = Ha((x)), 

S((x)) = Hr2{o)). 



and (A-15) becomes 



Ar„in = v7(r, + 7'„)(8fir)-'/2. 
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b. Goldstein [1955, 1956] 

For voltage sinusoidal modulation and sinusoidal correlation Goldstein [4, 5] obtains 

< = 4a^,^V^. (A-16) 

Goldstein's work corresponds to Ti = T2 = °K and 

a^ = kTn, 
where k is Boltzman's constant. Thus (A-16) becomes 

c. Bunkin and Karlov [1955] 

For square-wave modulation and sinusoidal correlation, Bunkin and Karlov obtain [17] 






^ ^ AT ^ 1 AT^ 



T,,+Te 2(ri„+n; 



1/2 

(A-17) 



If the precorrelator bandwidth AHi is much greater than the postcorrelator low pass filter band- 
width Afl, then 



Using the following conversions 



AHi + AH 

^ in ^ n-, 

AO^tt/t, 



for M< 1, (A-17) becomes 



ATrnin = ^{Ti-^Tn)(2BT)-'i'^ 



d. Kelly, Lyons, and Root [1958, 1963] 
For square-wave modulation and detection, (49) and (31) of Kelly, Lyons, and Root [12] becomes 



AT=2[n + T„]yj^-^^^- (A-18) 
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Equations (A— 11) and (A-14) are used to establish the conversion for B and r. Using the following 
conversions, 

T2 = T2 = T,, 

r(0) = 0, (for zero gain fluctuations) 
B = B, 



(A-18) becomes 



Ar„i„ = 2(7', + 7'„)(«T)^'/^ 



e. Knight |1962| 

For rectanjiular modulation and rectangular correlation. Knight [9] obtains 

S_ 2jra{(Tl-a-i,f /_a\ / ^\ 

N~ -ZVat [a + hJV a + bj- ^^ ^^' 



wlicrc the 2tt omitted in Knight's eq (4) is added. Tlie appropriate resolution criterion is SjN = 1. 
With this condition, then 

al{a + b) = (L 
T = 7r/r, 
a = B. 
and (A-19) becomes 



AT^,,= Tn[Brd{l-d)]-^l^ 



8. Appendix B. Theorems on the Fourier Transtorm 

Several theorems involving the Kourier transiorm of the lime function yj{t) and jzit) of zero 
mean are proved in this appendix: 

Lim = 0, (B-1) 
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for yxit) uncorrelated with yz(t), and 



Lim ^ y..WFV,^(i.') =0, (B-2) 



unless y^(0 = yi(0, 7^(0=7^(0, and v'=v\ and 



_ .. Y{0) 

^^\^1^~r~^ (B-3) 



T- 



and 



and 



and 



and 



and 



Um^^ = j^ W,{v)dv, (B-4) 



Yx{v + pvq) * Yziv-V-qvo) = Yjcz{v + pi^o + <7^o), (B-5) 



1 re/2 

^/,r = o a^{t)eJ'^''^fdt. (B-6) 

yj-e/2 



/ j, r — ^ ^i, m^i, r- m^ (t) 7) 



QiM = ^P,,„.Pt„„ (B-8) 



where Pi, r and /I i,p are defined in (18) and (12) respectively, and 0=^Pq^ is one period of the periodic 
signal p(t), and 



'2CaPu~a = c{t)p,{t) (B-10) 

a 

where 

1 fe/2 

c(Opi(0^o c{t)p,(t)dt. (B-11) 

o j-e/2 

8.1. Theorem 1 

Taking the inverse transform (,^) of the left-hand side of(B-l), then taking inverse transform of 
Y\{p) and Y2(^^) and interchanging the order of integration, one obtains (Lim T^ ^ is understood) 

T~' [" e^^'^''Yi(p)Yt(i^)dp = T-' f f y,(t')y2(t')dt'dt" f" eJ2^^<^-*'^^"M^. 
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The integration over v yields the Dirac delta function 8(t — t'-f t") {24\. 
Thus integrating over t', 

T ' f e^2vr.tY,(^)Y*(^)d^' = T-i p y,(t+t")y2(l")dt". 

//yi(t) and y2(t) are not correlated, then the integral on the right is finite, and when divided by T 
tends to zero as T-^oo. Thus, the left-hand side is zero, the Fourier transform is also zero and 
{B-1 ) is proved. 

8.2. Theorem 2 

The theorem is proved in a manner analogous to theorem 1. 

8.3. Theorem 3 
By definition, 

y = Lim 7^ y(l)dt. 

T-'-oo 1 j-T/2 

The function y(t) is truncated so that y(t) equals zero for |t| > T/2. Thus, 

1 r- 
y = Lim;j; y(t)dt. 

Taking the Fourier transform o/y(t) and interchanging the order of integration, 

T--* 1 J -00 J-x 

r/ie integration over time yields the 6-function of v. Subsequent integration over v yields (BS). 

8.4. Theorem 4 

Let Yxz(^) be the Fourier transform o/yx(t)yz(t) where yx(t) and yz(t) are random variables with 
zero mean: 



(.)=/_ 



Yx.(«^)= e-^2-'yx(t)yz(t)dt. 

J -co 

// yx(t) i*5 no? correlated to yz(t) or e"^^^"*, then the mean of the product is zero. Thus, Yxzli^) is 
finite and if divided by T, will approach zero a5 T-> 00. IfVzM ~- YxM- ^hen 

J - T/2 

where the limits of the integral are taken from — T 12 to T/2. This does not change the value of the 
integral because y{i) = for |t| ^ T/2. When y|(t) is rewritten as y| + [y|(t) — y|] and the above 
equation is divided by T {the limit T —^ ^ being understood), 

f— fT/2 rT/2 _ ] 

T--^Yxx(i^)-=T-Mv| e-'^^^^'dx-h e-^2.^^[y2(t)-y2]dt • 

I J-T/2 J-T/2 J 
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The integral in the second term is finite and the term goes to zero as T ~^^. The first term on the 
right-hand side goes to zero unless v = 0. For v = the term approaches y|. Thus, 

fyl forp = Oandx = z, 

r->oo ""^^ ' [0 otherwise. 

Using (10) yields {B-4). 

8.5. Theorem 5 

Using the definition of convolution, eq (5), 

Taking the Fourier transform and interchanging the order of integrations. 

The integration over v' yields 6(t — t'), then the result of integrating over l' is identified as 
Yxz(^+mz^o + qi^o). 

8.6. Theorem 6 

To obtain {B-6), both sides of {18) are multiplied by e^^^^^'o^ and then integrated over a period, 
using for Pi(t) the equivalent a?(t) as noted in {14). 

8.7. Theorem 7 

Using {B-6) and (12), 






Due to the orthogonality of different Fourier components, the integral has nonzero value only when 
s = r — m. This leads directly to {B-7). 

8.8. Theorems 8, 9, 10, and 11 

The proof of [BS), {B-9), and {B-10) are analogous to the proof of {B-7). Equations {13), 
through (21) are needed for various definitions. 

9. Appendix C. Radiometer Using an IF Amplifier 

A schematic for the IF amphfier is shown in figure C-1. Instead of the ampHfier shown in 
figure 1, now a mixer, local oscillator, and IF amplifier are used. The amplifier noise temperature 
T„ of figure 1 is now interpreted as the noise temperature of the mixer-IF amplifier as referred to 
the input of the mixer. Using the same notation convention as the body of the paper, 

yLo{t) — rj cos 2Trv\t, 

ysM{t) = r]ys{t) cos 27rvit, (C-1) 

yA{t) = hA{t)irysM{t), 
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where T) is the amphtude of the local oscillator of frequency V\. From (C-1), it is clear that an 
ideal balanced mixer is assumed. Taking the Fourier transform. 



Applying (6), 



Ysy,{v) = h'n\Ys(v+v^)+Ys{v-v^)l 
YAv) = HAv)Ysm{v). 

WAv) = GAv)Ws,i(v). 



(C-2) 
(C-3) 

(C-4) 
(C-5) 



9.1. Calculation of }/ 2 

The calculation of yT^ proceeds in a manner similar to section 3.5 except (C-2) and (C-3) are 
used in place of (29). Thus, instead of (46), now we have, 



Substituting (C-6) into (30) will ^ive terms of the type [see eq (47)], 



(C-6) 



YAv + mv(^ ± Vx)-k YAv + qV() ±V\) = Yj,Av + mi^o + qv^) ± 2vx). 



(C-7) 



Using an argument similar to that used to obtain (48), noting tbat the terms like (C-7) do not 
contribute, 



Lim rp ~oV /,^i,q^i, -q-mi^i 



(C-8) 



Equation (C-8) is identical to (48) except for the multiplicative constant \ if. Thus the remainder 
of the calculation is unchanged except that the gain functions X of (51) now contain \ yf. Thus, 



X = |T?^A/y,2(0)[26',(0)*G,(0)]i/2. 
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(C-9) 



9.2. Calculation of (7^2 

The calculation of cr/, 2 proceeds in a manner similar to section 3.5 except (C-4) and (C-5) are 
used in place of (37). Because Ws{i^-^Pi) = Ws(v—Pi}, (C-4) and (C-5) combined become, 

Thus the calculation is similar except W'^sii^) is replaced by ^ rf^JFsii^}- If (C-9) is used for the 
definition of X, then the statement of >7 2 and orF2 in (50) and (64) are otherwise unchanged. 

9.3. Comments 

The bandwidth B is defined in a similar manner whether the amphfier is at the signal frequency 
or an IF frequency. For an IF amplifying system, B is the bandwidth of the IF amphfier. If for 
some reason the signal were only in the signal channel of the IF amplifier (the noise remains in 
both channels), then Ti and Tz would be replaced by 2 7^1 and 2 7*2, but no other changes. Thus, 
the single channel IF radiometer would be equivalent (as far as sensitivity is concerned) to a similar 
double-channel IF radiometer except it would have twice the amplifier noise temperature. 
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